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Abstract. Recently, the semiclassical theory of the anomalous Hall effect 
induced by the Berry curvature in Bloch bands has been introduced. The theory 
operates only with gauge invariant concepts, that have a simple semiclassical 
interpretation and provides a clear distinction among various contributions to 
the Hall current. While the construction of such an approach to the anomalous 
Hall effect problem has been long sought, only the new semiclassical theory 
demonstrated the agreement with quantitative results of rigorous approaches 
based on the Green function techniques. The purpose of this work is to review 
the semiclassical approach including the early ideas and the recent achievements. 
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1. Introduction 

The anomalous Hall effect (AHE) is one of the most famous transport phenomena in 
magnetic materials. Unlike in paramagnets, the Hall resistance R xy of a magnetic film 
has two contributions. One is usual, it is proportional to the applied magnetic field 
H . The other one is anomalous, it is observable only in a ferromagnetic state. The 
anomalous contribution is often proportional to the magnetization rather than to the 
applied magnetic field 

Rxy = r H + r a M, (1) 
where M is the magnetization of the sample, ro and r a are constants that characterize 
the strength of the standard and the anomalous Hall resistivities respectively. The 
recent theoretical research demonstrated that the linear dependence on M in is 
not universal and the Hall resistivity can show resonance features as a function of 
variable parameters, including the magnetization. 

The practical interest in the AHE has continuously been driven by a difficulty 
to measure the carrier density in ferromagnets. Standard techniques based on 
measurements of the Hall conductivity are obstructed by the considerable anomalous 
contribution which can be much greater than the standard Hall conductivity. Recent 
advances in spintronics, especially the creation of new types of diluted magnetic 
semiconductors, revived the interest in the AHE as a useful tool to control spin- 
polarized currents and to characterize the magnetization. In addition, the recent 
theoretical interest has been fueled by the new interpretation of the anomalous 
Hall conductivity in terms of Berry phases and topological defects in the crystal 
band structure. Many theoretical constructions that usually had been considered 
of relevance mainly in high energy physics such as noncommuting coordinates and 
magnetic monopoles, became useful and even measurable in experiments on the AHE 

n la cam. 

Despite the long history and the considerable practical importance, the theory of 
the anomalous Hall effect has remained controversial. The first steps to explain the 
AHE in ferromagents were made more than 50 years ago. Since then many articles 
were published to correct previous mistakes and to suggest new explanations. Many 
of such efforts still were incomplete. While they resolved several pieces of the puzzle 
they also disregarded the others. Sometimes distinct quantitative predictions followed 
from applications of different methods to the same model. Such a controversy persists 
even at the present time. For example, there is a number of recent publications with 
contradictory quantitative predictions for the AHE in the Rashba 2D electron system 
HI C3 HI HI HOI [HI HH H3 HH ; although the issue has finally been resolved [H El • 

Many aspects of the AHE have been extensively reviewed in the literature. The 
detailed up to date discussion of experiments with diluted magnetic semiconductors 
and the comparison with existing numerical and theoretical predictions can be found 
in recent reviews [HI El El E] • The modern topological interpretation of the AHE 
and the Kubo formula in terms of a magnetic monopole in the momentum space 
was reviewed in Q15] . There are also much older introductions [201 [H] , concentrated 
on the side-jump effect in III-V semiconductors, although many concepts, discussed 
there, have been strongly revised in recent years. The recent review of these older 
work, made in the same spirit, can also be found in |22j together with the discussion 
of results on the spin Hall effect. 

Numerous efforts to design a rigorous semiclassical approach that would explain 
in simple terms all possible contributions to the anomalous Hall conductivity including 
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disorder effects, however, still lack a detailed comparison in a single work. The 
goal of the present review is to fill this gap and to discuss in more details the 
existing semiclassical theories, their advantages and limitations. We start with the 
earliest ideas introduced by Karplus, Luttinger and Smit and end with the most 
recent constructions that demonstrated the 1-1 agreement with the rigorous quantum 
mechanical techniques. 

The structure of this review is the following. In the rest of the introduction 
in section 1.1, we recall the basic text book information about the semiclassical 
approach to conductivity calculations taking as an example free electrons interacting 
with elastic scatterers. In section 1.2, I remind several commonly known facts about 
Bloch bands and explain how the introduction of the band structure complicates the 
creation of the semiclassical theory of the transport, especially in the application 
to the AHE. In section 2, I discuss the forces driving the AHE and review the 
earliest theories, including the quantum analog of the Boltzmann equation applied by 
Luttinger (section 2.1) and the introduction of noncommuting coordinates by Adams 
and Blount (section 2.2). In section 3, I proceed with more recent theories based 
on the gauge invariant formulation of the wave packet dynamics and its Berry phase 
interpretation. In sections 3.1 and 3.2, I will review the application of this approach 
to the anomalous Nernst effect, the intrinsic contribution contribution and the side- 
jump effect. Sections 4 and 5 are devoted to the rigorous semiclassical theory of 
the AHE, free of most limitations of previous approaches. In section 4, I discuss 
the rules that connect the scattering matrix with the classical concepts such as the 
scattering probability and the size of the coordinate shift at a scattering event. In 
section 5, I introduce the semiclassical Boltzmann equation and subsequently explain 
all important contributions to the Hall conductivity. There I will also discuss the 
strength of the AHE and comparisons with rigorous quantum mechanical approaches. 
Section 6 is the summary that discusses the present status of the theory and outlines 
possible future research directions. 

1.1. Semiclassical approach to conductivity calculations. 

Quantitative estimates of the dc AHE by standard techniques based on the evaluation 
of Green functions and their products invariably involve long complex calculations. It 
is hard to achieve transparent interpretations; therefore theories of the AHE normally 
focus on particular simple model Hamiltonians and ignore many-body interactions 
apart from mean-field exchange potentials that encode the magnetic order. Even with 
these simplifications, the AHE theory remains difficult to develop. 

One of the problems is the small magnitude of the AHE in comparison to 
the longitudinal conductivity. When considering the perturbative expansion of 
the conductivity in the weak disorder limit, the AHE contribution appears only 
in subdominant terms of higher powers in small parameters. Many standard 
approximations turn out to be no longer valid at these orders. Even a proper 
counting of relevant terms of a similar strength was a problem in many cases. Another 
difficulty is in the physical interpretation of conductivity contributions in the Kubo 
formula or in the Keldysh technique. Generally these rigorous quantum mechanical 
approaches operate with nongauge-invariant objects, such as off-diagonal elements of 
Green functions, of the density matrix or of the velocity operator, which only in the end 
are combined in the gauge invariant expression for the conductivity. Such calculations, 
while formally rigorous, hide the physical origin of elementary microscopic processes. 



CONTENTS 



4 



This complicates the analysis and the bookkeeping of the relevant contributions. 

The alternative approach is based on the classical Boltzmann equation applied to 
the electron transport [UJ (Ml 023 . It can be justified by the fact that in sufficiently 
clean materials one can look at the transport from the basis of wave packets rather 
than Bloch waves. In the dilute disorder limit, a wave packet is not destroyed during 
long time and behaves in many respects as a classical particle. One can trace the 
motion of the wave packet in external fields and describe it in terms that have a clear 
meaning in classical physics. 

In crystals, electrons cannot be considered as free particles because they strongly 
interact with a periodic crystal potential. As we will see in following sections this 
creates interesting ingredients in the wave packet dynamics but initially, we will 
describe the semiclassical theory free of these complications assuming that electrons 
do not interact with the crystal potential and with each other [23] . The impurity free 
Hamiltonian of such an electron system has plain wave eigenstates 

Mr,t) = j^e^-^ t , (2) 

where L is the size of the system, D is its spatial dimension and k = |k|. To construct 
a wave packet with a well defined average momentum k c , plain waves ([2]) should be 
superposed with the envelope function a(k), sharply peaked near the point k = k c so 
that / d D k\a(k)\ 2 k = k c , then the wave packet vector can be written in the coordinate 
representation as follows 

f d D k ( ( k 2 t \ 1 

*ke(r,<)=y ^a(k)esp|i(k.r- — Jj. (3) 

The normalization condition requires that 

(*kj*k c > = J d D rn c (r,t)* kc (r,t) = J d D k\a(k)\ 2 = 1 (4) 

and the index k c tells that the wave packet has this average momentum, namely, 
switching to the momentum representation one can find that 

k c = (* kc |k|* kc ) = J d D kj d D k' a (k)a*(k')<k'|k|k> = J d D k|a(k)| 2 k. (5) 
The velocity of the free wave packet center of mass can be derived as follows 



r|$ kc ) = d.{f<grj d D kfd D k'a(k)a* (k')e- jk ' r (re lkr ) e 1 ^-^ } 
Tt {jd D kJd D k'a(k)a*(k') (/ ^e-^ HJ^ kr )) e*^-<&} = 
|{/d D k/ci 1? kV(k / )5(k-k / )e i ^ 1 i| E [a(k)e- i ^]} = 



Tt {J d D k\a(k)\^} + | {/d°ka(k)' (^a(k))} = fc. 

In the external uniform electric field, the Hamiltonian operator is H — I r/2 + eE-f 
and the average of the momentum is changing with time 

k c = |(* kc |k|* kc } E = f t {j d D kf d D k'a(k)a*(k'){k'\e iJtt ke- iAt \k)} = 

(7) 

/d D k/d £, k'a(k)a*(k')(k'| [k, -ieE-rl |k) = -eE, 
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where k and f are respectively quantum mechanical momentum and coordinate 
operators. From (JT]) it follows that under the action of only the electric field the wave 
packet will accelerate indefinitely. This never happens in metals because of scatterings 
on impurities, that randomly change the direction of motion. It is impossible then 
to trace trajectories of all wave packets and the natural language to describe such a 
system is provided by the semiclassical Boltzmann equation. 

In classical physics the Boltzmann equation is the evolution equation for the 
particle distribution in the phase space. We will always assume in this work that the 
system is spatially uniform on scales much larger than the distance between scatterers, 
where the classical Boltzmann equation for scatterings on elastic impurities has the 
following form [H[25] 

^-^ = -£^(*-/w. (8) 

k' 

The rhs of ([5]) is called the collision term. For electrons that interact only with 
static impurities but non with each other the collision term is a linear functional of 
the distribution function. This linearity is not affected by the Pauli principle [27J. 
However, when many body interactions contribute to the collision term the Pauli 
principle leads to contributions proportional to /k(l — /k') e.t.c. We will not consider 
the latter case in our discussion. The scattering rate Wk,k' depends on details of 
the scattering potential and should be found separately. For a sufficiently smooth 
impurity, one can use wave packet equations to find the scattering cross-section 
by purely classical means but, in most of realistic applications, a smooth potential 
approximation does not hold for an impurity. Often the opposite limit of a (5-function 
type of a potential is considered as a reasonable assumption. This fact jeopardizes the 
applicability of the semiclassical approach but quantum mechanics provides a simple 
solution. There is the rule that connects the quantum mechanical scattering matrix 
with the classical scattering rate. This rule is called the golden rule of quantum 
mechanics. For a weak impurity potential in the lowest Born approximation it reads 

m 

Wk.k' = y|Vk,k'| 2 S(ek-e k '), (9) 

where is the kinetic energy of an electron with the momentum k and Vk,k' 
is the matrix element of the disorder potential between two states of an electron 
before and after the scattering. In what follows we will assume that % = 1. The 
potential of randomly placed impurities is V(r) = 5Zj^(r — Rj), where i enumerates 
impurities, R, are their random positions and v(r — R) is the potential of the single 
impurity with its center placed at Rj. One can show [27], that for such a disorder 
(|Mc,k' \ 2 )dis — ^|wk-k'| 2 ! where n is the impurity concentration and Vq is the Fourier 
transform of the single impurity potential at R; = 0. 

Together, the golden rule © and the classical Boltzmann equation ((5J) allow to 
perform the quantitative self-consistent calculation of the conductivity. Assume that 
the electric field is weak and look for the solution of the Boltzmann equation in the 
form 

/k = iW(e k ) + <?k, (10) 

where / eg (ek) is the equilibrium distribution, which depends only on the energy of a 
particle and does not contribute to the current and c/k is the correction linear in E. At 
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the steady state the term with the partial time derivative is zero and to linear order 
in E the correction to the distribution satisfies the time-independent equation 

d fegjtk ) 
<9e k 



eE ■ v k ^^ = - £ c^Gfr - (11) 



k' 

where in our case v k = k/m. Looking for the solution in the form <? k = g E • k, one 

can find a self-consistent result 

er tr E-k df eq {e k ) 

5k = — 5 , [12) 

m oe k 

where 

— = ^ w k,k'(l-cos(k,k')), (13) 

T tr k , 

r tr is called the transport life time. The electric current is given by the expression 

J = -e^ 3k v k . (14) 

k 

Using that at zero temperature Ql^tiSel = — 8{e F — e k ) one can arrive at the following 
expression for the conductivity along the electric field in the 2D electron system 

e 2 k F v F T tr , 1C , 
<Jxx = ; ■ (15) 

The important point is that the result (|15p is rigorous in the sense that 
when conductivity calculations are performed by formally exact quantum mechanical 
techniques, such as the summation of disorder averaged Feynman diagrams in the 
Kubo formula in the diffusive regime and disregarding higher order effects, such as 
the weak localization, one arrives at the same quantitative result. 



1.2. Difficulties with semiclassical approach in application to the AHE. 

The above discussion of free electrons demonstrates that it is possible to derive 
transport coefficients using the classical Boltzmann equation. The power of this 
approach is in its transparency. It operates only with concepts that have simple 
classical interpretations. There is no problem with the gauge invariance. 

Because of its simplicity this approach is ideal for introducing to the physics of 
the electron transport. Many concepts of solid state physics can be explained with 
a sufficient rigor without using complicated Green function techniques. Also, the 
semiclassical approach is needed to develop the scientific intuition about the model. 
Having done calculations for a simple problem one can be interested in further more 
complicated phenomena. It is always good to have a preliminary expectation about 
the final result. The semiclassical approach allows to make such an insight, while it is 
considerably harder with other techniques. 

Certainly, the semiclassical approach has limitations but it should not be 
considered as valid only in the classical limit. It makes rigorous estimates even for 
impurities that have no analog in classical physics, such as for delta-function potentials. 
Hence, its domain of validity is larger. However, in applications to the AHE, the 
semiclassical theory faced with a number of complications. Sometimes it has been 
speculated that the AHE is a purely quantum mechanical phenomenon that cannot 
be explained by classical means [55]. We will discuss later that it is not true but first, 
I explain the main arguments that created this scepticism. 
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Figure 1. The 2D electron system described by the Hamiltonian (1171 . 



In real crystals electrons are not free particles. They interact strongly with the 
periodic potential of the lattice. It has been well known that the lattice periodicity 
does not result in random scatterings. The Bloch theorem guarantees that eigenstates 
of the electron Hamiltonian in a perfect crystal have the form 

VWr, t) = ^e^-^VkCr), (16) 

where u n k(r) is a periodic in the elementary unit cell function. Due to this 
simplification, one can describe other interactions, such as scatterings on lattice 
imperfections with an effective Hamiltonian, not dealing directly with the lattice 
potential. The price for this is that the wave function has generally a nontrivial 
periodic part u n k(r) and the dispersion is no longer quadratic, et ^ k 2 /(2m) 
(usually it is possible to approximate it near important symmetry points by a 
quadratic dependence with a renormalized mass). The spectrum is also not everywhere 
continuous and splits into bands. 

The multiple band structure plays a very important role in the theory of the 
AHE. Because of it, operators of many observables are matrices in the band index 
space, which can have nonzero off-diagonal (inter-band) elements. One of the simplest 
models of this kind, shown in Fig [lj is the Rashba coupled 2D electron system with 
an additional out of plane Zeeman interaction. Its Hamiltonian reads 

H = — + \(k y (T x - k x a y ) - Aa z , (17) 

where &i are Pauli operators. This Hamiltonian describes two bands with different 
dispersions 

e ± = ^- T V(Afc) 2 +A2, (18) 
2m 



CONTENTS 



where +/— stand for major/minor band indexes. The velocity operator v = dH /dh 
has following components 



k x /m i\ \ „ / ky/ra A 
— iX k x /m J ' y I A ky/m 



(19) 



It is straightforward to check that neither v x nor {j y commute with the Hamiltonian. 
This means that in the Bloch basis of eigenstates of the Hamiltonian (fTT)) , the velocity 
operator still has nonzero off-diagonal elements, for example, 

if- ) , v+-* 0, „-+ ^ 0. (20) 

The semiclassical interpretation of diagonal (intra-band) matrix elements of the 
velocity operator is trivial. If one prepares a wave packet made of Bloch states of 
one band, the free motion of such a wave packet will be the corresponding diagonal 
velocity. The off-diagonal velocity matrix elements are more subtle. They do not 
affect the motion of a free wave packet. Only if a coherence among states of different 
bands is introduced due to some perturbation can their expectation values become 
nonzero. This mixing can be produced e.g. by an applied external electric held. Due 
to the Bloch vector dependence of the periodic part of the Bloch wave, the coordinate 
operator generally has nonzero inter-band matrix elements 

(u n k|r|ttnk) = {Unk\i~^-) ^ 0, (21) 

where n ^ n and (Mnk|w n 'k'} is understood as the integral of the product of periodic 
parts of Bloch states over the unit cell or, in the case of the Rashba Hamiltonian (fT7|) . 
it means the scalar product of spinor parts of Hamiltonian eigenstates. 

When the electric held is applied, the total Hamiltonian has a contribution 
He = eE • f , that not only accelerates wave packets, but also mixes states of different 
bands. Because of this, the expectation of inter-band parts of velocity operator 
components becomes nonzero. In turn, this means that in the applied electric held 
the instantaneous velocity of the wave packet is no longer a corresponding diagonal 
part of the velocity operator, but contains an extra (anomalous) component. The 
situation reminds the chiral anomaly in the quantum held theory [29j where the 
noncommutativity of the axial current operator with the Hamiltonian leads to effects, 
unexpected from the Hamiltonian symmetries. 

Similarly to the electric held, the impurity potential also mixes states of different 
bands. For example, the point-like impurity potential V(r) — Vod(r) in the chiral 
basis of the Rashba 2D electron gas has the following matrix form 

vo ( KK) («+!«£) 



^ = 75- k T k k ' , (22) 

with nonzero off-diagonal matrix elements (see [12 for explicit expressions). This 
means that the impurity role does not reduce to a simple instantaneous change of a 
direction of the particle motion. When a wave packet passes near such an impurity 
its wave function becomes distorted and the inter-band part of the velocity acquires a 
nonzero expectation due to the local band mixing. Thus in the vicinity of impurities 
wave packets move along unusual trajectories. 

The AHE was found to be related exactly to that type of microscopic processes. 
Hence the construction of the semiclassical theory of this effect faced with the problem 
of how to include the inter-band coherence into the purely classical description. 
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2. Early theories of the AHE 

The first theoretical proposal to relate the AHE and the spin-orbit interaction was 
made by Karplus and Luttinger [30]. They started from the fact that due to the 
relativistic corrections the effective Hamiltonian of an electron in a periodic lattice 
potential V(r) has an extra contribution due to the spin orbit interaction 

Hso = -4^2<HP x W). (23) 

This part of the Hamiltonian modifies Bloch wave functions of electrons, introducing 
specific Bloch vector dependence in their periodic parts itnk(r). When the electric 
field is applied the corresponding term in the Hamiltonian 

H E = eE-f (24) 

has nonzero matrix elements between states of different bands 

(u n k\H E \u n u) = ieE-(u„ k |^^) ^ 0, (25) 

ok 

where n ^ n. This band mixing ultimately leads to an unusual linear in the electric 
field contribution to the transverse velocity. 

The spin orbit coupling alone does not lead to the AHE because the anomalous 
transverse velocity, even if present, would have different signs in different degenerate 
bands unless the time reversal symmetry is broken. In ferromagnets this symmetry 
breaking appears spontaneously due to the exchange interaction, which is often 
approximated in theoretical models by the mean Zeeman-like field acting on electrons 
spins. In the simplest picture of the AHE this Zeeman field creates a population 
imbalance between bands with opposite signs of the anomalous velocity thus leading 
to the Hall current, proportional to the magnetization. This scenario is not always 
correct. For example, the recent research showed that often the strong AHE appears 
near points of the spectrum where the magnetization field lifts degeneracies in the band 
structure [19] . leading to resonance features at the Fermi energies in their vicinity. 

Karplus and Luttinger did not expect that their theory would be the final answer 
to the AHE puzzle. The reason was that the current they found was not gauge 
invariant, and thus could not describe the real observable. They pointed, however, 
that the magnitude of the found expression and its dependence on the impurity 
concentration were in a good agreement with results of experimental measurements, 
thus suggesting that their theory captured basic microscopic physics behind the AHE 
and the spin-orbit coupling and the magnetization must be the forces, driving the 
effect. 

Their theory was subsequently criticized by J. Smit |31[ I32j. who made the first 
effort to design the gauge invariant theory of the effect. Smit agreed with that the spin 
orbit coupling is responsible for the AHE but suggested a different mechanism. His 
approach was semiclassical in spirit. Smit's suggestion was to look at the evolution of 
a wave packet, as a semiclassical object and to design a Boltzmann-like equation to 
describe the evolution of the wave packet distribution function in the phase space. 
Today his name is usually associated with only one of the contributions to the 
Hall current, called the skew scattering, which was a new effect, not discussed by 
Luttinger at that time. However, in addition to the idea of the skew scattering, Smit 
discussed other contributions. Tracing the evolution of a wave packet he also found 
the anomalous velocity in the external field [32], though he did not think that this 
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velocity contributes to the Hall conductivity because, he thought, it is exactly canceled 
by another effect. 

Smit pointed that the anomalous velocity follows from the change of the 
polarization of wave packets by showing that the average of the coordinate of the 
wave packet contains generally an additional component A due to the Bloch vector 
dependence of the periodic part of the Bloch wave function 

r c (k,t) =v k t + A = —t+ (u k \i—u k ). (26) 

The second part A = (uk|«^Uk) depends on k. According to Smit, when the electric 
field accelerates the wave packet, the vector k changes and consequently this changes 
the polarization A. The wave packet becomes deformed. This evolution of the 
polarization leads to an additional charge transport in the transverse to the electric 
field direction. Smit pointed that A is not gauge invariant and rather its curl 

F = curl k A (27) 

should enter the final result. In the modern terminology A and F are called the Berry 
connection and the Berry curvature respectively. 

Smit's objection to the relevance of the anomalous velocity to the AHE 
conductivity also deserves a discussion. Smit pointed that in the DC limit wave packets 
cannot be constantly accelerated. While the electric field changes the polarization by 
accelerating wave packets, scatterings on impurities produce on average the exactly 
opposite change of k if the system reaches the steady state. Thus Smit concluded that 
coordinate shifts at scatterings should have exactly opposite effect on the wave packets 
polarization and thus on the AHE conductivity. The coordinate shifts at scatterings, 
first introduced by Smit [31] , indeed are the important part of the modern AHE theory. 
They were named "side-jumps" by Berger [33j[34l[35] who studied the effect in more 
details. 

Smit's work was the precursor of the modern semiclassical approach. He made 
the first effort to understand the anomalous Hall conductivity in classical terms such 
as corrections to the velocity of wave packets, coordinate shifts at scatterings (side- 
jumps) and asymmetric scatterings at an impurity potential (skew scatterings). All 
these ideas are currently incorporated in the theory, although his conclusion about the 
exact cancellation of the intrinsic and the side-jump contributions is not supported by 
rigorous calculations. 

There are two main reasons why his arguments fail. One is that the " polarization" 
A is not a good quantum number and, in fact, is not gauge invariant because it changes 
under an arbitrary momentum dependent change of the phase in the definition of Bloch 
states; therefore one cannot apply classical balance arguments to it. The second point, 
omitted by Smit, is that the side-jump can lead to the asymmetry of the distribution 
function even without an asymmetry in the collision term kernel in the Boltzmann 
equation. Such a distribution asymmetry is rather due to the change of the kinetic 
energy that particles experience after the side-jump in the presence of an electric field. 
The corresponding correction to the distribution function was named the anomalous 
distribution [53]. When coupled to the conventional part of the velocity de^/Ok the 
anomalous distribution leads to the Hall current. 
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2.1. Luttinger's rigorous theory: the quantum Boltzmann equation. 

In 1958 Luttinger published the detailed study of the AHE [33] based on the 
rigorous quantum mechanical approach that he had designed with Kohn in a previous 
publication [27] • Later this approach was generalized by Lyo and Holstein to the 
regime of ac external fields [37]. Luttinger's theory was correct but it did not find 
the general acceptance as a calculation tool and later many researchers have been 
looking for alternative techniques. The reason was that Luttinger's approach is very 
nontransparent. It involves many equations that self-consistently determine nongauge- 
invariant values. 

In this section I will try to make an explaining introduction to the Luttinger's 
paper. Rather than to follow directly his steps the goal here is to show that it is 
possible to explain the Luttinger's derivation with a simple schema and a different 
notation, according to which it is easy to classify AHE contributions in the Luttinger's 
approach and to make connections with the scmiclassical theory 

Luttinger starts with the evolution equation for the density matrix 
dp 
dt = 

where H is the Hamiltonian that includes both the disorder part and the electric field 



(28) 



H = H Q + V(r) + eE x x. (29) 
In the stationary state one should require that 

| = 0. ,30, 

If the solution of |28|) is found, the transverse electric current is given by the 
expression 

J y = -eTr [v y p] . (31) 

Hats mean that objects are matrices in the band index space. Since the velocity 
operator v y is diagonal in the momentum space, only the momentum-diagonal part of 
p is needed to calculate the current. However, v y can have off-diagonal elements in 
the band index space. From (|28[) . ([3H]) Luttinger derives the analog of the Boltzmann 
equation, which contains terms that depend only on the diagonal in the momentum 
space part of the density matrix. Schematically, it is useful to group terms that appear 
in his quantum Boltzmann equation as follows 

E[DT}{p eq ) - i[H ,p„eq} = Icol(pneq), (32) 

where Hq is the part of the Hamiltonian, independent of the electric field and of the 
disorder potential (Luttinger worked in the basis of Bloch states that diagonalize Hq). 
p eq is the equilibrium part of the density matrix (in the presence of the disorder but 
in the absence of the electric field) and p neq is the correction, linear in E. [DT] means 
the "driving term" which explicitly couples to the electric field. In the Bloch basis 
the driving term can be written as a series in powers of the disorder potential V, that 
starts at V°. 

[DT] {p eq ) = [DT] (0) + V 2 [DT] (2) + . . . . (33) 

The collision term for elastic scatterings on static impurities is linear in p neq and 
also can be written as the series that starts at V 2 . 

Icol(Pneq) = V 2 I { V{p neq ) + V 3 I { c 3 J(p neq ) + V 4 4% neq ) + .... (34) 
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This suggests to look for the solution for the nonequilibrium part of the density 
matrix in the form of a series in powers of V . Separating terms of the same order in 
V we find the chain of equations. The first equation allows to determine the largest 
term in the expansion of p neq and others allow to express higher corrections through 
the lower ones. From (|32|) . (|33[) and ([34)) and the linearity of I co i as a functional of 
Pneq it follows that this series begins at the term of order V~ 2 , i.e. 

Pneq — V Pneq > v Pneq * v Pneq T • • • < (OO) 

Simple power counting shows that to determine the correction of order V~ 2 it is 
enough to keep the driving term at zeroth order in V and the first term in the collision 
part of (|31|), i.e. 



E 



([DTr) =(/<2(pU a >)) , (36) 

\ / diaq \ / diaq 



where the index diag means that we take only the band-diagonal part of the expression. 

( — 2) 

Luttinger found that pheq is diagonal in the band index and does not contribute to 
the Hall current. It, however, makes the dominating contribution to the longitudinal 
current and is needed for further calculations. 

Next order contribution p^ 1 ^ satisfies the equation 

^]fc 1 )+4olfc 2) ) = 0, (37) 
where p^~ 2 ^ is already found by solving (1361) . 

It turns out that pneq , found from (|3"T|) . is still diagonal in band indexes and 
contains the antisymmetric contribution in the transverse to the electric field direction. 
It leads to the transverse conductivity that, like pneq , depends as 1/n on the impurity 
concentration n. 

At zeroth (next) order in the disorder strength, both inter-band and intra-band 
matrix elements become important. One can separate 4 distinct parts. 

Pneq = Pint + Psj + padist + Psk, (38) 

where first two terms are purely off-diagonal and the other two are diagonal in band 
indexes. These contributions to the density matrix satisfy following equations. 

({DT}^) -i[H ,p int ]=0. (39) 

V 7 of f -diag 

Note that terms in ()39[) do not depend on the impurity potential. 
The equation for p s j reads 

-i[HoJ. j ] = (l$(ptt ) j) ff „. • (40) 

v ' of f -diag 

p S j is purely off-diagonal in band indexes and appears because the collision term 

I^lipneq) has a nonzero off-diagonal part. 

The next contribution is diagonal and follows from the compensation between the 
higher order driving term and the collision part. 

E[DT]^=I^{p adist ). (41) 

Finally, there is a contribution due to the compensation between two collision 
terms 

€(Ps k )+C(pL?) = 0- (42) 



E 
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All four parts of the nonequilibrium density matrix in (1381) contribute to the 
Hall current via (|3"Tj) . Interestingly, since they are formally of zeroth power in V and 
the velocity operator does not depend on V the resulting conductivity due to these 
contributions for a Gaussian correlated disorder becomes independent of the impurity 
concentration. 



2.2. Noncommuting coordinates. 



Luttinger's theory is both complete and well justified quantum mechanically. However, 
his approach features the same problems as other rigorous quantum mechanical 
techniques in applications to the AHE. It does not explain what is happening in 
simple terms. What can be concluded from his work is that off-diagonal elements of 
the density matrix and of the velocity operator play the important role. Separately, 
they are not gauge invariant and only their product produced finally a gauge invariant 
current. Because of these complications, Luttinger did not discuss the semiclassical 
meaning of the derived contributions. 

To resolve this issue, in 1959 Adams and Blount made an effort to create the 
semiclassical theory based on the introduction of noncommuting coordinates [38] . 

The straightforward semiclassical approach, based on the preparation of a wave 
packet from states of the same band, may fail when the electric field is applied. Since 
the electric field mixes states of different bands, a part of the initially free wave packet 
starts fast oscillations with frequencies to ~ e n ^ — e fi .k in comparison to the rest of 
it. Such a wave packet does not satisfy the basic criteria of being a classical object 
because it would be composed of parts with strongly different oscillation frequencies. 

The resolution of this problem was first suggested by Adams and Blount |38j . 
The off-diagonal part of the Hamiltonian due to the electric field can be considered as 
a periodic field that modifies the Bloch wave functions. Thus one can choose another 
Bloch basis, in which the term with the electric field has no inter-band matrix elements. 
To linear order in E = E x x the periodic part of such modified Bloch states reads 



4k) = Kk} + ieE : 



\ ' ( - "nk|%T x k ) 



Eftk 



Kk), 



then at t = one can prepare a new wave packet 

d D k 



<k» 



L D / 2 



a(k)exp{i(k-r)}\u' nk ), 



(43) 



(44) 



which in the external electric field does not split right away into differently oscillating 
components. It is now instructive to calculate the velocity of the wave packet in the 
transverse to the electric field direction. Let the electric field points along x-direction 
and H = Hq + eE x x is the full Hamiltonian (we do not consider impurities yet). Then 
the wave packet evolves according to 



< kc (r,t)= e -^< kc (r). 
The transverse velocity is calculated as follows 



Kk) 



de„k i A p Tp ( I du„ k dunk\ _ I du„ k i du„ k \\ _ de„ k 
dk v ^ lcrjx \\ dk v I dk x I \ dkx I dk y I ) ~ dk y 



(45) 



(46) 
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where F™ is called the Berry curvature of the Bloch band with index n and 
is its z-component. 

The first term in the last line of (f46|) is just the usual velocity that equals the 
diagonal part of the velocity operator, while the rest of the expression is called the 
anomalous velocity. This anomalous contribution is ultimately responsible for the 
intrinsic AHE. 

There is sometimes a misunderstanding about how the noncommuting coordinates 
appear in the theory. Sometimes it is simply stated that this happens to coordinate 
operators after the unitary transformation that switches to the basis of Bloch states. 
This, of course, is not true because the unitary transformation alone cannot make 
commuting operators noncommuting. The reason is more subtle. In the standard 
Bloch basis with the periodic part \u n u) that diagonalizes Hq the coordinate operator 
has the form 

r = z|^ + A(k)+X, (48) 

where A(k) is diagonal in band indexes and its elements are the Berry connections of 
Bloch bands and X is purely off-diagonal in band indexes. If we want to work with a 
modified Bloch band with periodic parts of Bloch states given in (|43p we should switch 
to that new modified basis, i.e. to make the additional rotation which transforms \u n k) 
into \u' nk ). In that basis Hq is no longer pure diagonal and contains the off-diagonal 
component that cancels with eE • X. Thus in the modified basis the full Hamiltonian 
has the form 

H = H + eE-(i-^+A(k)), (49) 

where now Hq is, in fact, new operator, that has the same matrix form in the modified 
basis (|4"5|) as the old free Hamiltonian Hq in the original Bloch basis |it„k)- 

In the basis (|43|) the Hamiltonian ([49]) is by construction diagonal in band indexes 
including terms linear in the electric field. Thus, when working with modified Bloch 
states, it is useful to regroup the terms in the original Hamiltonian so that in the 
modified basis it looks like we still have the unperturbed band-diagonal part H but 
instead of the usual coordinate operator the electric field couples only to its projection 
on the subspace of a given band. The same observation holds if we consider the average 
of any operator over the state of the wave packet (]44[) . prepared from the modified 
Bloch states 

(* nkc (r,t)|f|*„ kc (r,t))«^'(k-k c )+A n (k c ), (50) 

where we used the fact that the inter-band component of the coordinate operator has a 
zero expectation value and higher order terms in the electric field can be disregarded. 

Summarizing, if we consider the evolution of the wave packet made of modified 
Bloch states, mathematically, instead of working with band mixing by the electric 
field, we can assume that we still deal with the original free diagonal Hamiltonian Hq 
but coordinate operators should be substituted by their projected versions 

f-fc = «^+A(k). (51) 
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The anomalous velocity appears now as the result of the noncommutativity of 
such modified coordinate operators. Thus for the electric field along the x-axes the 
transverse velocity contains a component 



v 



y = -i[y c ,eE x x c ] = -ieE x 



1 dk„ + A y , i gk + A 3 



- P E (d^L-dA* 
crjx \ dk x dk y 



(52) 



Adams and Blount pointed that their approach can also be applied to the 
scattering problem on the impurity potential if the latter is sufficiently smooth so 
that a moving in its field wave packet feels only a weak gradient. One can truncate 
the off-diagonal part of the impurity potential V(r) — > V(r c ) so that the wave packet 
acquires an additional velocity 

v? = -i[? i c,H]^j- (53) 

The comprehensive discussion of the noncommuting coordinates in Bloch bands can 
be found in [39]. The theory of Adams and Blount was applied by Chazalviel |40j . 
Nozieres, Lewiner [41], Berger [33 [35 EH] and by Lyo and Holstein [42] to the AHE 
in III-V n-type semiconductors. Because of the degeneracy of electronic bands in 
their model, they had to extend the theory to the case, which today is called "bands 
with a nonabelian Berry curvature" . They started from the standard 8-band model 
with a spin orbit interaction. Due to the hybridization the conducting bands are also 
influenced by the spin orbit coupling. Similarly to the above discussion, Chazalviel and 
Nozieres projected coordinate operators to the space of only two conducting bands. 
The problem, however, is that the latter have the same energy dispersions and thus 
the subsequent separation of conducting bands is impossible and one has to keep the 
off-diagonal elements of the coordinate operator in this subspace. Thus the Berry 
connection and the Berry curvature become 2x2 matrices. 

d 

r c = i— + As x k, (54) 
ok 

F = 2As, (55) 

where s is the pseudospin operator acting in the space of the conducting bands index. 

This nonabelian case, however, is trivial because the Berry curvature does not 
depend on the Bloch vector. If one considers a 2D electron system, only the conserved 
out of plane component of the pseudospin enters the Hamiltonian and the problem 
reduces to two separate bands with the effective impurity potential 

V( r ) -> V(r) + As- (k x W(r)) . (56) 

The second term in (|56l) is functionally similar to the spin orbit coupling due to 
relativistic corrections ([2"3"|) but with a different strength A, which can be considered 
as a renormalized spin orbit coupling constant. Formally, corrections to the impurity 
potential and to the coordinate operator ([B"4")l and ([5"6"]) should be included in the full 
Hamiltonian of the 2D electron gas, together with the Rashba coupling. However, 
their effects are usually weak because of the weakness of the parameter A in ([5"6"]) , that 
makes the corresponding Hall conductivity also small 40J. In contrast, the Hall effect 
due to the Rashba coupling is nonperturbative. For example, the estimates of the 
intrinsic contribution for the Rashba coupled electron gas show the strength close to 
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the universal value e 2 /2h in some range of parameters [5]; so the effect of the Rashba 
term on the AHE is expected to dominate, although at large Fermi energies there can 
be a crossover between effects of two types of spin orbit couplings. 

The weakness of the Hall effect due to the coupling A in (|56]1 can be understood 
in terms of topological defects in the band structure. Nonzero A is induced by a weak 
mixing of p-type orbitals to states of the conducting bands; hence it is suppressed, 
approximately inversely proportionally to the cubic power of the large gap between 
conducting and valence bands [40]. In contrast, the Rashba coupling, even when small, 
creates the topological defect centered directly inside the conducting bands. At Fermi 
energies close to this "resonance" point the AHE is nonperturbative in the spin orbit 
coupling and hence is very strong. 

Despite a partial success of the Adams's and Blount's approach, their semiclassical 
theory has strong limitations. One is the difficulty to apply it to a short range impurity 
potential. Such a potential destroys the wave packet and thus cannot be treated 
in a weak gradient approximation. The second problem is that this approach still 
operates with not strictly classical concepts such as noncommuting coordinates. This 
complicates the interpretation of other objects in the semiclassical theory such as 
the distribution function. It also complicates the derivation of the skew scattering 
contribution. Thus in all publications following this approach, the important part of 
the skew scattering contribution was missing. That part is parametrically similar to 
the side-jump contribution and cannot be disregarded. For example, the Chazalviel's 
and Nozieres's et al conclusion that the total impurity concentration independent Hall 
conductivity in their model is the same in the magnitude but has opposite sign from 
the intrinsic contribution is wrong because of this omission. 



3. The Berry phase theory of the AHE. 

Sundaram and Niu [43j designed a very powerful and unifying framework to study 
the wave packet kinetics. Their approach is based on the derivation of the effective 
Lagrangian of a wave packet moving in weak fields. The idea of the approach is 
that the time-dependent Schrodinger equation for a wave packet is realized from the 
variational principle with the Lagrangian given by 

L = <*kji^ - #1*0. (57) 

The time dependence of the wave packet implicitly is contained in the time dependence 
of its average momentum k c and coordinate r c = {^k c |r l^kc) and possibly in other 
explicitly time-dependent parameters in the system. This allows to rewrite the time 
derivative in (|57p in terms of r c and k c 

, • /, \ .* / i du i . . . du . . .du, ,„_. 
L = k c r c - e(k c ,r c ) + k c (w| — ) + r c (u|^> + (58) 

where \u) is the unit cell periodic part of the Bloch wave. 

The first two terms in the Lagrangian Q58J) are the same as in a typical Lagrangian 
of a classical particle with the classical mechanical Hamiltonian H c i — e(k C) r c ). The 
other terms are geometric, in the sense that their contribution to the action depends 
on the trajectory in the phase space but not on the rate of the motion along this 
trajectory. This is the feature of the Berry phase in quantum mechanics and effects 
of the last three terms in (|58|) can be called Berry phase effects. 
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The Lagrangian formulation provides a fully gauge invariant approach to the 
study of the wave packet dynamics. It is now easy to find the equations of motion 
in an arbitrary potential with a weak gradient in a magnetic field B by varying the 
action over the trajectory. For the motion in a band with a dispersion ej. the result 
reads 

^ c = dk" ~ * 

(59) 

k c = -eE - W(r c ) - er c x B, 

where F is the Berry curvature of the Bloch band. Wave packet equations show 
that the Berry curvature can be considered as an unusual magnetic field acting in 
the momentum space. However, unlike the magnetic field in electrodynamics, which 
is a pure curl, the Berry curvature originates from a source. In electrodynamics, a 
magnetic field with such properties would originate from magnetic monopoles. Their 
analogs in Bloch bands are the points of exact crossings of band dispersion curves 

Bang. 

The wave packet equations can be generalized when the motion in degenerate 
bands is considered. If states of degenerate bands mix coherently the evolution 
becomes more complicated. The state of the wave packet should then be considered 
having finite amplitudes in both bands [44] [45] 

|<& kc ) = /" dkaC^) [*7i|*i> + »ftl*2>] , (60) 

where (i — 1,2) are basis functions of the wave packets in each band. The 
coefficients r\i enter the effective Lagrangian and should be considered as independent 
variables. Their dynamics was found to be according to equations 

dm ( ■ , . duj \ 

,_=^-k c ( Ui |t— )j Vj , (61) 

where are matrix elements of the Hamiltonian. 



3.1. The anomalous Nernst effect 

The recent triumph of the wave packet approach was the first semiclassical explanation 
and the quantitative theory of the intrinsic anomalous Nernst effect (ANE), which 
is the AHE driven not by an external electric field but rather by the temperature 
gradient. The theory of the ANE [46] is based on the previous understanding of the 
intrinsic angular momentum |43[ I47j of a wave packet and also on the observation |48j 
that the wave packet equations ([59]) lead to the specific expression for the phase space 
volume. When considering the continuous limit it reads 

J2 - / [dk] (1 + eB • F) . (62) 
k J 

The physically measurable transport current is defined by 

j = J-VxM(r), (63) 

where J is the microscopic current and M(r) is the magnetization density. The latter 
can be found from the grand canonical potential in the magnetic field 

F= --^log(l + e- /5(eM ^ ) ) =-- /[dk](l + eB-F)log(l + e- /5(£M ^ ) ) , (64) 
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where (3 = 1/kgT and tu — e(k) — m(k)-B is the electron energy in the magnetic field 
coupled to the magnetic moment of a wave packet m(k) = — i(e/2)(V k it| x [H(k) — 
ek]|V k u). 

The magnetization is the magnetic field derivative of the thermodynamic potential 



"<->--(£) 



(65) 



The expression for the intrinsic Hall current, induced by a weak temperature gradient 
originates from the phase space volume correction in (|64[) . Substituting (|65l64p into 
(63) one can find [46] 



Jint 



VT 

-e x 

T 



J [dk]F [(e M - M )/ (e k )A :B Tlog (l + e^^^)) 



(66) 



5. 2. Wavepacket theory of intrinsic contribution and side-jump effect. 

The intrinsic contribution to the AHE is a straightforward consequence of the 
anomalous velocity term in wave packet equations (|59p . Under the action of the 
electric field, all the electrons in the Fermi sea will shift in the transverse direction, 
leading to the intrinsic Hall current [1] 

J'«* = -e 2 Ex /[dk]/ k F, (67) 

where /k is the electron distribution function in the given band and F is the Berry 
curvature. 

Sinitsyn et al demonstrated |7 that wave packet equations can also be successfully 
applied to the problem of a scattering on an impurity if the latter has a sufficiently 
smooth potential, thus providing the fully gauge invariant theory of the side-jump 
effect. Integrating ([59)) over the time interval during which a wave packet "feels" the 
impurity potential and assuming that this potential is sufficiently weak, one can find 
a coordinate shift at a scattering event 

K- k = Fx(k'-k), (68) 

where k' and k are center of mass momentums of the wave packet respectively after 
and before the scattering. This definition of the anomalous coordinate shift at a 
scattering on an impurity (the side-jump) is different from some expressions suggested 
in the early theories. For example Berger [211 HI] assumed that Sy a ~ k x . The 
Bcrgcr's definition follows from the identification of the coordinate shift with the Berry 
connection. Since the Berry connection is not gauge invariant, the old definition does 
not have a direct semiclassical meaning. In contrast, the expression (f6"5)) as well as 
more general expressions (|74p , (|75[) from next section are gauge invariant and depend 
both on the in-going k and the out-going k' Bloch vectors. This difference becomes 
important when constructing the rigorous theory of the effect when the scattering on 
an impurity is not isotropic. 

There are two main rather distinct effects due to the anomalous shift. One is the 
side-jump accumulation. After the averaging over many scatterings, side jumps do 
not cancel and lead to the velocity renormalization by a correction 

^'>(k) = 5> k ^<, !k . (69) 

k' 
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The second effect is that when a scattering takes place in the presence of an external 
electric field, there is a change in the potential energy upon a scattering given by 

At/ k >, k = eE • K<,k- (70) 

This change of energy ultimately influences the Hall conductivity and should be 
properly included in the semiclassical Boltzmann equation. 

According to [7], the side-jump related conductivities depend only on parameters 
taken near the Fermi surface. This is in contrast to the intrinsic contribution that 
depends on the integral of the Berry curvature over the whole Fermi sea. This may 
be one of the reasons why the comparison with experiments showed a good agreement 
with the intrinsic contribution calculations that disregarded impurity effects, except 
keeping a finite life time of quasiparticles |4J. Although the intrinsic and the side-jump 
contributions to the transverse conductivity do not depend on impurity concentrations 
and in this sense are parametrically similar, still if the Berry curvature is weak near 
the Fermi surface but strong deep inside the Fermi sea, the intrinsic contribution can 
dominate. 

The ultimate example of this kind is the quantum anomalous Hall effect that 
appears when the Fermi level is placed inside the gap in the bulk spectrum. In such 
insulators the gapless excitations are forbidden (except near the sample edges), so the 
side-jump and the skew scattering effects do not contribute to the conductivity but 
if the band has a nonzero Berry curvature there is a quantized intrinsic contribution 
to the Hall current. This happens e.g. in 2D Dirac bands, related to the graphene 
system [49l [50] . 



4. Scattering rules. 

The "philosophy" of the semiclassical approach is to operate only with classical 
concepts, however, using several rules that connect some of them with purely quantum 
mechanical ones in order to achieve a quantitatively rigorous result. The expression 
for the Berry curvature is one of such rules, namely it relates the anomalous velocity 
to Bloch wave functions. The scattering is described in quantum mechanics by the 
scattering matrix, which has no analog in classical physics; therefore one should use 
the rules that connect the scattering matrix to the classical microscopic effects. 

One such a scattering rule is widely known, and for its importance it is named the 
golden rule of quantum mechanics. Lets introduce a combined band-momentum index 
I = (n, k). The golden rule relates the scattering rate low from the state V into the 
state I in the continuous spectrum with the corresponding element of the scattering 
T- matrix [26]: 

low = 2tt\T w \ 2 5{e v - e,). (71) 

The T-matrix is defined as 

Tiv = <Z|V^'>, (72) 

where V is the impurity potential operator and \ifii) is the eigenstate of the full 
Hamiltonian H — Hq + V that satisfies the Lippman-Schwinger equation 

m = \l) + ir—^tyl). (73) 

ei - Hq + ir) 

Scattering rates low cannot include all the possible information encoded in the 
scattering matrix. This is obvious because entries of the T-matrix are complex 
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numbers and entries of the matrix of scattering rates are real. In the golden rule only 
the absolute value of the T-matrix elements are represented. Thus the semiclassical 
approach, which uses only the golden rule as relating the classical and the quantum 
descriptions of the scattering, should generally fail. This is indeed the case in the 
AHE. The golden rule does not contain the information about the side-jump effect at 
a scattering event. This fact forced authors of [51] to search for the gauge invariant 
expression for the side jump that would connect it to the scattering matrix. 

Such an expression indeed can be derived. In the lowest Born approximation it 
has a particularly simple form, 

d d 

5rm = (up\i—-u v ) - (ui\i—ui) - T)^^arg{Vi^i), (74) 
ok' ok 

where arg[a] is the phase of the complex number a and 

- d d 

Dk ' k -ov + av 

This type of expressions has been first found even before the work [51] . however, 
beyond the AHE theory. Belinicher et al derived it to apply in the photovoltaic effect 
[52] . They showed that when electrons absorb a polarized light they make shifts 
(|74p . where V would be responsible for the electron-photon interaction. Such shifts 
of the form (|74p finally contributed to the photo-induced conductivity in their model. 
Unfortunately the work [52] had been unnoticed by the modern Hall effect community 
until the expression ([74]) was independently rederived in [51] and applied to the AHE- 
problem. 

The expression for the side jump (|74|) is gauge invariant. Interestingly, it restores 
the information, lost in the golden rule. Unlike the golden rule that in the lowest 
Born approximation depends on the absolute value of the scattering potential, the 
coordinate shift expression depends on its phase but does not depend on its absolute 
value. Thus the expression (|74[) can be considered as complimentary to the golden 
rule. 

If the impurity potential is spin independent, then the side-jump does not depend 
explicitly on the type of the impurity potential and can be expressed in terms of initial 
and final states only [5TJ [53] : 

d d 

Srm = (uii\i——ui*) - (ui\i—ui) - D k /. k arg[(u//|u ; )]. (75) 
ok' OK 

At a weak scattering angle it reduces to equation (|68p . derived in previous section 
[4"9"] . The derivation of the golden rule ([7T]) and the expression for the coordinate shift 
(|75p can be done by considering a scattering of a wave packet. Imagine a state, that 
initially coincides with the Bloch state ^ (r, t) under the influence of a weak potential 
of an impurity V(r). The solution of the time-dependent Schrodinger equation can be 
written in terms of the eigenvectors ipi> (r, t) of the unperturbed Hamiltonian 

^ out (r,*)=S;C,/(t)#(r,t). (76) 
v 

Consider the wave packet, that was initially gathered around the state Then in 
the lowest order in the strength of the potential the perturbation theory leads to the 
following expression for time-dependent coefficients CV(i) (see Eq. 19.9 in [54] ) 

C v {t) = -iV vl f e^'-^'dt' + 5 n . (77) 
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The higher order terms can be incorporated into the above formula by merely 

substituting the T-matrix instead of the disorder matrix elements (see Eq. 19.10 
in [51]). 

Cu(t) = -%T V i [ e^'-^'dt' + 6m. (78) 



From this solution one can show that for I' ^ I, Cy (t) change with time according 
to the law 

C ^^=2«\T l , l \ 2 5(e v -e l ). (79) 

The coefficient |CV(i)| 2 has the meaning of the probability of the electron to be in the 
state V, from which immediately the golden rule ([71) follows. 

The coordinate shift expression was derived in a similar fashion. One can prepare 
a wave packet, approaching the impurity at the point Vi mp — according to the law 
r(i) = vo;i at t — > — oo. The scattering on an impurity generally destroys the wave 
packet, however having the scattering matrix one can write a formal expression for 
the average coordinate of the state of the wave packet after the scattering. According 
to |51j it can be written in the form 

r+oo= [ dr(*°"*(r,+ W ))*r* OTt (r,+c^) = ]T |C P (i ^ +c<))| 2 (v 0/ ,t + <5r r; ) • (80) 

Since \Cy{t — ► +oo)| 2 has the classical meaning of the scattering probability into the 
state I', the expression in the parentheses is reasonable to interpret as the coordinate of 
the particle if it is scattered in the state I'. Then the term Voj't simply tells that after 
the scattering the particle moves with the new velocity Voz' and thus the expression 
Srin in ([50)1 can be interpreted as the coordinate shift at a scattering event. This 
identification leads finally to the gauge invariant expression (|T4[) . 

The side-jump is a weak effect and thus its expression in the lowest Born 
approximation is sufficient for the further discussion. However the lowest Born 
approximation in the golden rule is inappropriate for the theory of the AHE. For 
a weak disorder one can use the expression of the T-matrix in terms of the Born series 
in powers of disorder potential matrix elements 

T,r«*fr+E VU " V Z ■ +-... (81) 

One can consider only several first terms in this series in order to capture the basic 
microscopic processes. Substituting Eq. ([8lj) into Eq. ([7T|) one will arrive at the 
expansion 

w u ,= U $+w$+u,$ + .~, (82) 



where 



( 2 ) - "WIU..I2 



IV 



2ir(\Vu>\ 2 ) dls 6(6 l -e l ,), (83) 



and so on. The skew scattering contribution to the Hall effect follows from the 
antisymmetric part of the scattering rate 

4^^^. (85) 
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Since u>if/ is symmetric, the leading contribution to u>ffi appears at order V 3 , at which 
the antisymmetric part of the scattering rate is particularly simple [36 ; 55 

colf, a) = -(2^) 2 £<5(e, - e vl )lm{V w V VV 'V vn ) dis 8{ei - e v ), ( 86 ) 
i" 

with the superscript 3a meaning that this is the antisymmetric part of the scattering 
rate calculated at order V 3 . Usually, properties of the skew scattering were inferred 
only from this lowest order antisymmetric part of ojw . Thus it is customarily assumed 
that w,,, is proportional to the impurity concentration. This, however, holds only 
in the lowest nonzero order, i.e. for w,,, . In the next order the antisymmetric 
scattering is proportional to the product of four disorder vertexes. For a Gaussian 
correlated potential (V • V ■ V ■ V)di S ~ (V ■ V)di S {V ■ V)dis ~ n 2 , where n is the 
impurity concentration. Thus the higher order contribution behaves as u>if> ~ n 2 , 
unlike wjj' ~ n - This means that the higher order term is different parametrically 
and has a rather distinct microscopic origin; therefore it should not be disregarded. 
Moreover, the contribution to the conductivity, arising from oj\f, , is parametrically 
similar to the side-jump related contributions. Hence, to derive all important effects 
one should calculate the symmetric part of low in the lowest Born approximation and 
then calculate its antisymmetric part, including next two orders in V. 

ujiv w w w + lj u , + u u , . (87) 

5. Mechanisms of the AHE in the semiclassical Boltzmann equation. 

Eqs. (|71[) and (|75p contain the quantum mechanical information necessary to achieve 
quantitatively rigorous predictions working in the framework of the semiclassical 
Boltzmann equation. This approach takes into account both the change of the 
direction and the coordinate shift during a scattering in a homogeneous crystal in 
the presence of a driving electric field E. Keeping only terms up to the linear order 
in the electric field the Boltzmann equation reads [53] 



dfi „ dfo(ei) ^ 
__eE. Vo| = _^ WJJ/ 



ft - Jv 5 eE • br w 

aei 



i' 

where expressions for ujui and 5ru> were derived in previous section, vq; is the usual 
velocity 

v 0( = dQ/dk (89) 

and fo(ei) is the equilibrium distribution. 

The Boltzmann- type equation ([55]) has the standard form, as (JSJ for free electrons, 
except the allowed inter-band transitions and the coordinate shift effect which is taken 
into account in the last term in the collision integral on the rhs of Eq. (|88|) . The 
derivation of this term has been explained in [3 |49l [53]. It follows from the fact 
that under the scattering in the electric field from the state /' into the state /, the 
side-jump is associated with the change of the potential energy AUw — eE • 8rw, 
which has to be compensated by the change of the kinetic energy. Thus the 
conservation of energy requires that e;/ — e; = eE • Srw . This in turn means that 
the equilibrium distribution does not annihilate the collision term anymore because 
fo(ei) - fo(ei>) ~ -^f^eE • dr w , which is exactly the last term in ([88]). The validity 
of these arguments was confirmed in numerical simulations ■ 
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The total distribution function ft in the steady state (dft/dt = 0) can be written 

as 

ft = fo(ei) + gt + gf + gf + gf ls \ (90) 

where gf, gf 1 , gf 2 and g adist are non-equilibrium corrections to the distribution 
function of linear order in the electric field (the label adist stands for the anomalous 
distribution). They solve self-consistent time-independent equations: [S3] 

-eE.v 0i ^ = -£^>( 5 ?- 5 ?,), (91) 

E 4 3a) (gt - g!> ) + E u w gf) = o, (92) 

I' v 

E 4< a) (9? - g s v ) + E 4< } igf -gf) = o (93) 



and 

E w «' (V Ist - gf lst - ^^ eE • ) = °- ( 94 ) 

One can deduce the dependence of the distribution corrections on the impurity 
concentration by noticing that ~ n, then from (p?Tj) follows that gf ~ n _1 . Then 

from (|52"]1 and the fact that ~ n follows that <7? a ~ then from w^ a ' ) ~ n 2 

and (J93J follows that gf a - n° and from (J94J) follows g; odlst ~ n°. A detailed solution 
of equations (|9H94|) in a special model was demonstrated in [53J . 

The first correction gf is symmetric and others gf dlst , gf 1 and gf 2 are asymmetric 
in the transverse to the electric field direction. When coupled to the usual part 
of the velocity the latter lead to 3 separately gauge invariant contributions to the 
conductivity. For the electric field along the x-axes these read 

< lst = -e^2(9t dlst /E x )(v ol ) y ~ n°, (95) 
l 

° s v f=-eY,(9t l /E x ){vo l ) y ~n- 1 , (96) 
i 

v S y? = -eJ2(9f/E x )(voi)y~n°. (97) 
i 

The last two contributions can be called the skew scattering conductivities because 
they originate from the asymmetric part of the collision term kernel low ■ It is, however, 
reasonable to distinguish between the two because of their different dependence on the 
impurity concentration. The first skew scattering conductivity (|96[) is the conventional 
skew scattering, that has been discussed by many authors 36, 55 . The second one 
(|97p was generally discarded, although it is parametrically the same as the side-jump 
conductivity. The explicit quantitative estimates of a^ 2 so far exist only for the 
massive 2D Dirac band [53J . Because of the lack of the proper terminology we will call 
(|97p the intrinsic skew scattering because this conductivity, similarly to the intrinsic 
contribution, is independent of the impurity concentration n. 
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The asymmetric corrections to the distribution do not exhaust all mechanisms 
of the AHE. We already discussed that between scatterings under the action of the 
electric field, wave packets move with an extra velocity 

vf = eExF;. (98) 

This velocity is linear in the electric field, therefore we did not consider its effect 
on the distribution function. However, when coupled to the equilibrium part of the 
distribution it produces a finite Hall current. 

f (ei)F, ~ n°, (99) 



Finally, the accumulation of coordinate shifts after many scatterings can be 
described semiclassically on average as an additional velocity contribution 

vf*=5^w,/i«r,, /> (100) 



n°. (101) 



Thus the total transverse conductivity can be written as the sum of five 
contributions: 

total int i adist , sj , ski , sk2 /ino\ 

a yx = < 7 yx+ a yx + a yx+ a yx + a yx ■ ( 1U ^) 



5.1. Semiclassical versus fully quantum mechanical techniques. 

Besides the semiclassical theories, the research on diluted magnetic semiconductors 
also stimulated the theoretical interest in other quantitative approaches. Due to the 
relative simplicity of several important models, such as the Rashba coupled 2D electron 
system, a number of publications appeared recently with rigorous quantum mechanical 
calculations by Kubo and Kubo-Streda formulas [H HOI H31 HH1 H31 1^1 [SZ1 [SHI ES] and by 
a variety of the quantum Boltzmann equation and the Keldysh techniques [81 151 |6"UI IBT] . 
Sinitsyn et al [51[ 153] demonstrated the 1-1 correspondence between semiclassical 
contributions to the AHE and the summation of relevant subseries of Feynman 
diagrams in the Kubo-Streda formula [62] ■ Similar agreement was established with 
the Luttinger's theory [36] . The results are summarized in Table [TJ 

According to 53J, the classification of contributions in the Kubo formula is not 
merely by a separation of diagrams into the disorder free part and the vertex correction 
but rather by the parts of the velocity matrices in chiral basis that stay inside the trace 
of the Kubo formula. Thus the intrinsic contribution appears from the summation 
of all diagrams with only off-diagonal parts of the velocity vertexes in Bloch basis. 
From this point of view the skew scattering is the most "classical" because it is due 
to the summation of all diagrams with only diagonal parts of velocity operators and 
the difference between conventional and intrinsic skew scatterings is due to different 
types of disorder vertexes involved. The conventional skew scattering is due to the 
vertex correction that involved correlators of three or more disorder vertexes while 
the intrinsic skew scattering is due to only Gaussian disorder correlations. The 
side jump and the anomalous distribution effects are related to Feynman diagrams 
that contain one off-diagonal (inter-band) and one diagonal (intra-band) parts of the 
velocity operator. This reflects the fact that although the side-jump itself is related to 
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Table 1. Mechanisms of the AHE in the semiclassical Boltzmann equation (SBE), 
the Luttinger's quantum Boltzmann equation (QBE) and the Kubo-Streda IG2I 
formula (KSF). Or . and v° d , stand for diagonal and off-diagonal parts of the 

x / y x j y 

velocity operator in the Bloch state basis, n is the impurity concentration, Tr is 
the summation over all states in the momentum space and over all band indexes, 
G R ,G A are respectively retarded and advanced Green functions of the electron 
system and other symbols have the same meaning as in the bulk of the text. 



mechanism 
of the AHE 


strength 
of CTyx 


SBE 


QBE 


KSF 


intrinsic 


O(n ) 


-eTr(/o«g) 


-eTr(p tnt v° d ) 


II | I ,int 
yx ~r a yx 


side-jump 
accumulation 


O(n ) 


~-eTv(g s v s y J ) 


-eTr( Ps3 v° d ) 


£ Tr (v° y d G R v$G A ) 


anomalous 
distribution 


O(n ) 


-eTr(g adist v 0y ) 


-eTr(p aMst v d ) 


fi Tr (v d G R v° d G A ) 


conventional 
skew scattering 


Oin- 1 ) 


-eTv(g al v 0y ) 

K, k ' - ») 


-eTr( / 3(- 1 )^) 


f^Tr {v*G*viG A ) 
(nongaussian vertex) 


intrinsic 

skew scattering 


O(n ) 


-eTr(g a2 v 0y ) 
( < k , - n 2 ) 


-eTr(p sk v d ) 


^(v d G*v d x G A ) 
(gaussian vertex) 



the anomalous velocity it contributes to the final current only after the electric field 
distorts the distribution function by a simple acceleration, i.e. by the coupling to the 
usual velocity in the Boltzmann equation or coordinate shifts create an anomalous 
distribution that again contributes to the final current via the coupling to the usual 
velocity. 

5.2. Terminology in the AHE theory. 

After the above discussion of all the effects leading to the AHE it seems appropriate to 
reconcile some of the differences in the terminology in recent publications. Generally 
it is stated that there are three basic microscopic effects, leading to the Hall current: 
the intrinsic contribution, the side-jump and the skew scattering. Table Q] shows that 
this classification is too restrictive. There are, in fact, five separately gauge invariant 
contributions, each having rather distinct origins, ft is possible to regroup them 
into three, because the side-jump accumulation and the anomalous distribution both 
originate from coordinate shifts at scattering events; similarly, the conventional and 
the intrinsic skew scatterings both appear from the asymmetry in the collision term 
kernel in the semiclassical Boltzmann equation. 

However, such a simplified classification into only three parts was one of the 
reasons for a confusion. For example, it is customarily stated that the skew scattering 
always leads to the conductivity that depends as 1/n on the impurity concentration. 
Table f shows that this is not true. The intrinsic skew scattering conductivity is 
independent of the impurity concentration and is parametrically very similar to the 
side-jump related contributions. Unjustified claims resulted in the omission of the 
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intrinsic skew scattering in almost all discussions of the semiclassical approach to 
the AHE. Such an omission has been repeated in the recent efforts to design the 
spin Hall effect theory by the analogy with previous AHE results [35]. It is also 
useful to distinguish the side-jump accumulation and the anomalous distribution 
effects. The side-jump accumulation is a rather direct consequence of coordinate 
shifts while the derivation of the anomalous distribution requires several extra steps 
in the semiclassical theory, which were unnoticed in a few former publications. 

Another confusing example from the recent terminology is the statement that 
the vertex correction, coming from a Gaussian correlated disorder in the Kubo 
formula is due to the side-jump effect only. Comparing the vertex correction with the 
semiclassical expression for the side-jump conductivity, discussed by Berger, Nozieres 
and others, one would find a discrepancy because the intrinsic skew scattering is also 
captured by the vertex correction and it was not considered by the older semiclassical 
theories that concentrated only on the side-jump effect. 

The classification into the "intrinsic" and "extrinsic" contributions was also 
understood quite arbitrarily by many authors. In the present review we coined the 
word "intrinsic" for the single special contribution, which is due to unusual trajectories 
of wave packets in the external electric field rather than any other mechanism that 
involves scatterings on impurities. This definition of the intrinsic contribution is 
justified by the fact that it follows only from the crystal band structure. Respectively, 
all other contributions can be called extrinsic. Sometimes, in other publications [19] 
any conductivity contribution, independent of the impurity concentration n is called 
intrinsic. This seems not a good choice of a terminology because the side-jump and the 
intrinsic skew scattering effects satisfy this definition but originate from scatterings on 
impurities. While corresponding conductivities are independent of n they can depend 
on other disorder parameters, for example, impurities with different typical ranges of 
scattering angles can result in Hall conductivities different by a numerical factor of 
order unity. 

There are examples where the "intrinsic" is associated with any effect induced 
by the Berry curvature in Bloch bands, while "extrinsic" would be due to relativistic 
corrections to the impurity potential, possibly renormalized in conducting bands by 
the crystal potential, as discussed by Berger, Chazalviel and others [33j [40l [41] . For 
example, according to this terminology all effects produced in the 2D electron system 
by the Rashba spin orbit coupling are intrinsic and all other effects such as due to the 
spin orbit part of the impurity potential are extrinsic. In such a case all the effects 
discussed in this review would be called intrinsic. Such a terminology also seems 
somewhat misleading, because effects due to the disorder spin orbit coupling also can 
be described and classified in the same way as here when working with the relativistic 
Dirac equation for electrons or, in the case of semiconductors, this corresponds to 
working with the 8-band model without projecting all operators to the conducting 
2-band system. Then one can work with a disorder potential free of the spin-orbit 
coupling and contributions to the AHE then can be derived in the same way as in this 
review [401157] , 

Finally, there is a notion of the " Berry phase" contribution to the AHE lj . Often 
it is identified with the intrinsic contribution as defined in this review. However, origins 
of all disorder related effects can also be traced to the Berry phase or to a nontrivial 
topology of Bloch bands. For example, the side jump expression and the antisymmetric 
part of the scattering rate can be expressed via a topological Pancharatnam phase [51J. 
Therefore it is possible to speak about the Berry phase mechanism of the AHE, which 
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includes all the physics discussed here but to apply this terminology to a particular 
contribution can be misleading. 

6. Summary 

The modern semiclassical theory rigorously takes into account all known important 
contributions in the model of electrons in Bloch bands interacting with static 
impurities. Predictions of this theory were verified with rigorous quantum mechanical 
techniques. However, so far the semiclassical theory of the AHE was built to deal with 
electrons that do not interact with each other. 

The state of the art is currently at the stage similar to where the theory of 
electrons in metals was before Landau introduced the Fermi liquid hypothesis. The 
Fermi liquid theory was originally semiclassical and allowed to derive many important 
properties of the electron state prior to systematic diagrammatic calculations. 

The important problem now is the effect of the nonzero Berry curvature on many- 
body interactions and only recently new publications appeared that addressed it. 
Thus, Haldane [63] proposed that the Berry phase can be considered as the property 
of quasiparticles living near the Fermi surface. In [64] Shi et al showed the robustness of 
several results of the wave packet theory against e-e interactions. Shindou and Balents 
in [65) studied the problem in more details by deriving the quantum kinetic equation. 
The interesting finding was that the Berry curvature now acts as a pseudo-magnetic 
field in the extended (k, u>) space. Also, Shi and Niu 66J considered interacting wave 
packets in bands with the Berry curvature and found the attracting force that can 
induce the instability in the p-channel and thus can facilitate the unconventional 
superconductivity. More generally, there is a strong similarity between properties of 
systems with the AHE and the p-type superconductivity and superfluidity in 3 He- A 
[67l IB"8l IB"9l [70l I7T] . The effective action governing low-energy physics of the p x + ip y 
superfluid contains a topological contribution, leading to the Hall effect, similar to the 
intrinsic contribution to the AHE in metals. It should be interesting to explore the 
diffusion equation for Bogoliubov's quasiparticles on this background. For example, 
one can expect to find effects similar to the side-jump and the skew scattering. 

Another poorly understood problem is the physics near the edges in systems with 
the AHE. It is possible that the side-jump type scatterings from the edge lead to 
the edge current. Analogous phenomena are known in geometrical optics where the 
reflected beam is shifted from the incident point on the mirror [731 [731 E2 ESI EH EH 
[75]. There can be complications with this analogy because a reflection of the beam 
generally changes its intrinsic angular momentum [78] . 

Recently the optically induced AHE attracted some interest both theoretically 
[El [79] [80] and experimentally 81, 82 . The conventional AHE needs a magnetization 
in order to break the time-reversal symmetry. This requirement can be avoided if 
other interactions, such as with a polarized light are introduced. An observation of 
such an optically induced Hall current would allow to explore the physics discussed 
above avoiding many difficulties in the interpretation of the standard AHE because 
currently ferromagnetic samples are very dirty. So far experimental results on this 
topic are controversial, with one group reported the observation of the effect [81] 
and another reported that the effect was not observed in GaAs at least up to the 
measurement uncertainty [53] with results in the agreement with the classical Hanle 
effect. The rigorous theory of the optically induced AHE is also missing, although the 
semiclassical theory allows to make an insight. For example, one can expect to find an 
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intrinsic-like contribution due to distorted trajectories of electrons in simultaneously 
applied DC and circularly polarized AC electric fields. The absorption of a photon 
from a polarized light can induce side-jump-like shifts etc. 

One more recent experiment demonstrated that the anomalous Hall conductivity 
can be measurable even in paramagnetic materials when electron spins are polarized 
by an applied external magnetic field [83] . The theoretical model of this effect must 
inevitably deal with a strong conventional Hall effect and the strict separation into 
conventional and anomalous Hall effects may not work anymore, for example due to 
the phase space volume correction following from wave packet equations of motion. 
The theory of the AHE in this regime should be upgraded. 

It is worth mentioning that new types of the AHE have been recently proposed for 
an experimental verification [73 [76l HH [85] and the semiclassical theory can be used 
to describe these effects too. There are also suggestions of alternative mechanisms of 
the AHE and the spin Hall effect based on the possibility of a spin-dependent force 
[5r?I [571 155] . This idea so far is based mainly on the Drude- model type of arguments 
and its verification by rigorous quantum mechanical and numerical techniques still has 
not been developed. 

Finally, the extra numerical and ab-initio study of the AHE is needed. So far all 
existing research of this kind concentrated only on the intrinsic contribution, treating 
discrepancies with experiments by introducing a finite life-time of quasiparticles [3ll89|. 
The rigorous numerical study of the anomalous Hall conductivity in models with 
a realistic disorder is still missing. In contrast, currently there is a number of 
successful publications on the related spin Hall effect (90J, [91] that demonstrated 
a good agreement with existing theoretical results and allowed to extend them to 
the analytically complicated and experimentally more realistic strong disorder cases 
[921 E3 Ell ISS [96l [971 [98]. Such numerical studies were also valuable to understand 
the problem of the spin accumulation near the edges. Similar research should be 
very valuable in applications to the AHE. Simple models such as the Rashba coupled 
2D electron system with an out-of-plane magnetization should be under the reach of 
modern numerical algorithms. 
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